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This paper presents an analysis technique for determining upper and 
lower bounds on the performance of a burst-trapping error control pro- 
cedure. The analysis is valid for random, burst, or compound channels 
provided that a block interleaving degree I ^ 1 can be found for the channel 
such that error patterns in blocks spaced I blocks apart occur approximately 
independently. Good agreement is achieved between the' theoretical per- 
formance of codes on telephone channels and the performance obtained by 
computer simulation. 

I. INTRODUCTION 

A burst-trapping error correction procedure for error control on 
compound channels such as the telephone channel has been described 
by S. Y. Tong. 1 An evaluation of the burst-trapping procedure by 
computer simulation of its performance on recorded telephone channel 
error data is presented in the companion paper. 2 An analysis technique 
for determining the performance of the burst-trapping procedure is 
presented here. The analysis technique permits determination of upper 
and lower bounds on the probability of block error. It is valid for random, 
burst, or compound channels provided that a block interleaving degree 
/ ^ 1 can be found for the channel such that error patterns in blocks 
spaced I blocks apart occur (approximately) independently. 

Expressions for upper and lower bounds on the probability of block 
error for codes of rate \ and f are given in Table I. For higher rate codes, 
numerical determination of the stationary probabilities which yield 
the bounds seems preferable to determination of expressions for the 
bounds. The performance of a rate § code, the (39, 26) shortened BCH 
code, has been computed using the expressions for probability of block 
error of Table I for two sets of recorded telephone channel error data. 
For both sets of data the computed performance agrees well with the 
performance obtained by computer simulation. 
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Table I — Bounds on Block Error Probability fob Rate 1/2 and 

Rate 2/3 Codes 

RATE J Codes 
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II. ENCODING AND BURST-TRAPPING DECODING PROCEDURES 

The encoding procedure and burst-trapping decoding technique have 
been described in detail elsewhere. 1 The code is a rate (6 — l)/6 recurrent 
code 3 whose parity check matrix, A, is constructed from the parity check 
matrix, H, of an (to, k) linear systematic block code [where k/n = 
(b — l)/b] and the (to — k) X (n — k) identity matrix I. The truncated 
parity check matrix, A N , for rate § codes is given in Fig. 1. The con- 
straint length of the code is N = [(& — 1)/ + 1]to where I ^ 1 is a block 
interleaving constant. 

Although encoding is specified by the A matrix, it is useful to interpret 
encoding as the interleaved encoding of I subcodes in the following way. 
Blocks 0, I, 21, 3(, • ■ • form the first subcode. The fc-tuple of information 
bits, Iu , of the [il]th. block (i = 0, 1, 2, • • •) is encoded into an n-tuple 
M n which is the concatenation of I it and a parity (n — fc)-tuple Q it : 



M it = /,-, ||Q„. 

I it can be represented as the concatenation of (b 
segments: 



J _ ri II j-2 II . . . II 7&-1 
* i t — * i t 1 1 L i t 1 1 1 1 ■* i t ■ 

The parity (n — /c)-tuple Q i{ is 

Qn = Pit + I{i-\)t + I{i-2)t + • • • + I(i-b+ni 



(1) 
1) equal length 

(2) 



(3) 



where P i( is the parity (n — /c)-tuple obtained by encoding I u with 
the block code parity check matrix H. Thus each encoded subcode block 
is a block code word whose parity portion is modified by the addition 
of an information segment from each of the previous (6—1) subcode 
blocks. Further the jth information segment of the [U]th block is added 
to the information portion of the [(i + j)t]t\i block (j = 1,2, ■ • ■ ,b — 1). 
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Fig. 1 — A n matrix of a rate 2/3 code. 
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A given block of a subcode is decoded by the burst-trapping procedure 
in one of three ways. Successive blocks of a subcode (blocks 0, 1,21, • • • ) 
are decoded by Random Error Decoding (RED) as long as the decoder 
decides that t or fewer errors occur in each block. When the decoder 
decides that more than t errors occur in a block, for example block 
[it], then that block is decoded by Burst-Trapping Decoding (BTD). 
The next (b - 1) subcode blocks (blocks [(i + 1)4 [(i + 2)4 ■■■ , 
[(i + b - l)t]} are decoded by Blind Faith Decoding (BFD). Detailed 
descriptions of RED, BTD, and BFD are given below. Successive sub- 
code blocks (blocks [{i + b)t\, \{i + b + 1)4 • • • } are decoded by RED 
until the decoder again decides that a block has more than t errors. 

Let Mf t , I? e , Q? { be the received n-tuple at the [i(\th block, the 
received information fc-tuple portion, and the received parity (n — k)- 
tuple portion respectively. Let Pf t be the parity (n — /c)-tuple obtained 
by encoding I? c with the block code parity check matrix H. Let l it be 
the decoded information /c-tuple of the [U]th block regardless of how the 
decoding is accomplished. 

The primary method of decoding is RED. RED is attempted at the 
[#]th block if none of blocks [(i - b + 1)4 [(t - b + 2)4 • • • ,[(»'- 
1)4 are decoded by BTD. RED is the removal of the effects of previous 
block information segments from the parity portion and the decoding 
of the block by bounded distance decoding of t or fewer errors. That is, to 
decode the [M]th. block the decoder computes 

PU - Qh + l\i-i» + Jcf-»« + • • • + JK»+«i (4) 

and decodes If t \\ P' it as a block code word. Following standard terminol- 
ogy for recurrent codes we will say that /}*-«< > *<*-«< > ' * ' > *(7-»+i>< 
are "fed back" to the [iflth block. 

The secondary method of decoding is BTD. BTD is effected when a 
RED attempt indicates that the block code word is detected in error 
but cannot be corrected by bounded distance decoding. If this occurs 
in block [U\ then blocks [(i + 1)4 [(*' + 2)4 •••,[(» + 6 - 1)4 are 
necessarily decoded by BFD (to be described subsequently) and decod- 
ing of block [U] is delayed until block [(i + b — 1)£] has been decoded. 
The information segments of block [ii\ are obtained from the corre- 
sponding parity portion of blocks [(i + 1)4 [(* + 2)4 • • ■ » [(* + 
b — 1)4 as follows 

*« = Qii + l)t T " (i' + l)/ I -La-Dt I -L(i-2U I " * * "T -»(i-6 + 2)/ | (5) 
*« = ^(i + 2)/ + ^(i + 2)/ + -l(i + l)/ T" -t(i-l)/ T " " " T" J(i-6 + S)/ ! (6) 
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I 3 il = Qti + »ii + Pfi+»)t + I(i + 2)t + I(i + l)t + I(i-l)l 

+ ... +I\:i b+iU ; (7) 

lit = Qfi + b-\)l 4" -Pfi + 6-lX 4" I(i + b-2)l "T I(,' + b-8)< 4" * * "4" i« + i)< ■ (8) 

In equation (6), for example, we will say that (Qf j+2) * + Pfi+nt) an d 
/(*«>< are "fed forward" and P (i . iu , •■■ , I b ( 7i b+au are "fed back" 
to the [U]th block. 

The third method of decoding is BFD. BFD is effected at the [U]th 
block if one of blocks [(i - b + 1)1], [(i -6 4- 2)4 ••-,[(»- 1)4 has 
been decoded by BTD. BFD of [#]th block is the use of the received 
information /c-tuple as the decoded information /e-tuple. That is 

I it m I* . (9) 

III. COMMUNICATION CHANNEL 

We consider channels where each block word (binary n-tuple) is 
subjected to the component- wise modulo-two addition of an error 
pattern (binary n-tuple) during transmission. We define a partition of 
the set of 2 n possible error patterns as follows. 

Co — is the set of one element — the pattern with no errors. 

C a - — is the set of channel error patterns which are identical to nonzero 

code words. 
C c — is the set of nonzero channel error patterns which are correctable 

by RED (bounded distance decoding of t or fewer errors). 
C u — is the set of channel error patterns which are both uncorrectable 

by RED and undetectable (but are not nonzero code words). 
C d — is the set of channel error patterns which are detectable (but not 

correctable by RED). 

We consider channels where error patterns separated by I or more 
blocks occur independently. On burst or compound channels, proper 
design of the burst-trapping procedure requires sufficient block inter- 
leaving that the requirement of independent error patterns is approxi- 
mately met. The nth power 4 of the binary symmetric channel with 
any I ^ 1 is included in the class of channels under consideration. 

Let p , Pw , Vc , P* > Pd be the probability that the channel error 
pattern for a particular block is in set C , C„ , C c , C u , C d respectively. 
Let a be the probability that the information bits of a received word are 
unaffected by the channel error pattern given that the channel error 
pattern is either correctable, detectable, or undetectable. For the class 
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of communication channels under consideration the probability that 
the error patterns in successive blocks of a subcode are in sets C, and 
d is PiPi , i, je{0, w, c, u, d\. 

IV. ANALYSIS OF PERFORMANCE 

Since each subcode is independently and identically decoded it suffices 
to analyze the decoding of one subcode. The decoding of successive 
blocks of a subcode can be described by a set of states and a set of 
transition probabilities between states which form a Markov Chain. 
The set of states is partitioned into a set of "normal" states and a set of 
"anomalous" states. Given the block code parameters, the sequence of 
states is determined by the error pattern sequence. The block code 
parameters which are required are d m , the minimum distance of the 
block code, and t, the amount of bounded distance random-error-cor- 
rection done in RED. 

Before enumerating the states we give a general description of the 
two sets of states. In all following discussion we assume the decoding 
of a single subcode. A block is decoded in a normal state only if all the 
blocks containing information segments fed back to that block are 
correctly decoded. Thus a block decoded in a normal state can only 
be in error due to its own error pattern or due to errors which are fed 
forward. A block is decoded in an anomalous state if one or more of 
the blocks containing information segments to be fed back to that 
block are incorrectly decoded. Once a block is decoded in an anomalous 
state the decoder is affected by fed back errors in addition to the channel 
error pattern sequence. Successive blocks are assumed to be decoded 
in anomalous states until a run of V blocks with no channel errors 
occurs. V is termed the recovery space for a subcode. V is determined 
in terms of d m and t (independent of details of the specific code) so 
that if no channel errors occur in a period exceeding (V — 1) blocks, 
then no errors are fed back to subsequent blocks regardless of the 
previous channel history. The first block after a run of V blocks free 
of channel errors is decoded in a normal state. A block decoded in an 
anomalous state can be in error due to its own error pattern, due to 
errors which are fed forward, and/or errors which are fed back. Error 
propagation which is due to fed back errors is bounded and can occur 
only in anomalous states. 

Each normal state is numbered by means of the following notation: 

The first digit represents the method of decoding 
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1— Random Error Decoding (RED), 

2— Burst-Trapping Decoding (BTD), 

3, 4, • • • , (b + 1)— Blind Faith Decoding (BFD). 

States whose first digit is 1 or 2 have a two digit number. States whose 
first digit is ;, 3 ^ j g (6 + 1), have a j digit number. The last digit 
is 0(1) to designate that the block is correctly (incorrectly) decoded. 
Interior digits for blocks with 3 ^ j ^ (6 + 1) digit numbers designate 
the decoding history of the (j - 2) previous blocks. An interior digit is 
0(1) to designate that its corresponding block is correctly (incorrectly) 
decoded. For rate (6 - l)/b codes there are 2 6+1 normal states. Normal 
state diagrams illustrating this notation are given in Figs. 2 and 3 
for rate } and f codes respectively. Throughout this paper a rate f code 
will be used in a running example. The analysis technique is easily used 
for any rate (b — l)/b code. 

Since the error pattern sequence determines the state sequence the 
transition probabilities between normal states can be expressed in terms 
of p , Pv, , Pc , Pu , Pa and a. Transition probabilities between normal 
states are given in Figs. 2 and 3 and Appendix A for rate \ and § codes. 
For simplicity of notation in writing the one step transition probabilities, 
p iti , from state i to state j the normal states are renumbered as in 

Figs. 2 and 3. 

The normal state diagrams represent the possible state transitions 
until an undetected error occurs (state 11) or a BTD and its associated 
(6 - 1) BFD's occur [states (6 + 1)00 • • • to (6 + 1)11 • • • 1]. Transi- 
tions from state (6 + 1)00 • • ■ or state (b + 1)10 • • • (states marked 
with in Figs. 2 and 3) are the same as those from state 10 since the 




x x 

Fig. 2— Normal state diagram for rate 1/2 codes. 
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(pu + Pw) 




Fig. 3' — Normal state diagram for rate 2/3 codes. 

last (6—1) blocks containing information block segments to be fed 
back are correctly decoded. Transitions from other normal states 
(marked with X in Figs. 2 and 3) are to anomalous states since one or 
more of the last (6—1) blocks containing information block segments 
to be fed back are incorrectly decoded. The enumeration of the anom- 
alous states is presented later. 



V. AN OPTIMISTIC ESTIMATE OF THE PROBABILITY OF BLOCK ERROR 

If we assume, for purposes of obtaining an optimistic estimate of 
performance, the presence of a Genie Decoder which corrects all fed 
back information, then transitions from normal states marked X (in 
Figs. 2 and 3) are the same as those from state 10 or the two states 
marked (f>. Note that with the Genie Decoder blocks can still be in- 
correctly decoded and incorrect information can be fed forward but 
only correct information can be fed back. The assumption of a Genie 
Decoder in effect eliminates error propagation (which arises from 
erroneous feedback) and the only possible states are normal states. 
The entire transition probabilitj' matrix [p,-,,] is then known from the 
normal state diagram (for example, Figs. 2 and 3) where p,,, is the one 
step transition probability from renumbered state i to renumbered 
state j. The Markov Chain is regular and therefore the stationary 
probabilities exist. 5 The stationary probability 7r, of being in renumbered 
state i, i = 1, 2, • • • , 2 +1 , can be obtained from 



2 » +1 
i = \ 



(10) 
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and 

E n = i. ( n > 

The probability of block error with a Genie Decoder, P ea , is the 
probability of being in a normal state whose identifying number has a 
last digit of 1. For rate f codes, 

P,a = *"s + ^6 + 1*8 + ^O + Tig + 7T 13 + 7T 14 + JT 16 . (12) 

Expressions for P., for rate \ and f codes are given in Table I. For 
higher rate codes numerical determination of the stationary probabilities 
which yield P eo seems preferable to determination of an expression 
for P e „ . 

5.1 Upper Bound on the Probability of Block Error 

An upper bound on the probability of block error can be obtained by 
considering worst case error propagation. The set of anomalous states 
is used to represent the decoder when error propagation can occur. 
Assume that one or more of blocks [(i — 6 + 2)4 to [it] is incorrectly 
decoded and that block [U] has channel errors but all succeeding blocks 
[(i + 1)4 [(i + 2)4, • • • have no channel errors. Error propagation 
is the effect that, although no channel errors occur beyond block [U], 
decoding errors in block [it] and/or previous blocks may, through feed- 
back, cause decoding errors to occur in some succeeding blocks [{i + \)l], 

[(» + 2)4, 

The amount of error propagation, a, is the number of blocks after 
the last block with channel errors that have decoding errors. The length 
of error propagation, X, is the number of blocks after the last block 
with channel errors up to and including the last block with decoding 
errors. Upper bounds A and A on a and X respectively are derived in 
Appendix B and given in Table II for codes with cl m ^ 2£ + 2. In Table 
II the notation \_xj means the greatest integer less than or equal to x. 

Another required quantity, the recovery space, is the minimum 
number of blocks after the last block with channel errors that must 
be free of channel errors to guarantee that error propagation ceases 
and to guarantee that the decoder has completed any BTD decoding 
and associated BFD decodings resulting from error propagation. 
Specifically, the recovery space, V, is defined as the minimum number 
of consecutive blocks free of channel errors (blocks [(i + 1)4, 
[(i + 2)4, • • • , [(i + V)(] required to guarantee that the decoder will 
return to one of normal states ! 10, 11, 20, 21) at block [(i + V + 1)1] 
regardless of the channel error sequence prior to block [(i + 1)1]. The 
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Table II — Recovery Space and Bounds on the Amount and 
Length of Error Propagation for Codes with d m > 21 + 2 



b = 2 



6 = 3, t = l 



a = l 
A = l 

7 = 2 



A = 2 
A = 2 

V = 5 



d„ > b 



b > 4, t = 1 



A = 36 - 8 

.4 = 6 

j/ = 46 - 7 



rf« ^ 6 



A = 36 - 8 

A - (6 - 1) + Ij^ + lj 

F = 46 - 7 



d m > bt 



bt 



b = 3, I = 2 



A = 36 - 6 

.4 = 6 

I' = 46 _ 7 



d™ ^ 3/ 



A = (26 - 3) 

I (6 - 2)t - 2 I 
L rf m - 2< J 

/I = (6 - 1) 

7 = (36 - 4) 

- 2)t - 2 



+ 



I (6 - 2)1 
L d m - 



2/ 



<L > 3/ 



A = 36 - 6 



A = (6 - 1) 



1(6-91-2 

^ L d m - 21 ^ 

V = 46 - 7 



recovery space for codes with d m ^ 2£ + 2 is derived in Appendix B 
and is given in Table II. Values for A, A, and V are presented in Table 
III for typical values of b, t, and d m . 

The decoder can be represented by (V + 1) anomalous states, as in 
Fig. 4, when it is not in a normal state. State A exists whenever channel 
errors occur in a block. States A x , A 2 , • • • , A v represent successive 
blocks free of channel errors. In Fig. 4 state X represents any normal 
state marked X in Figs. 2 and 3 from which the transition from normal 
state to anomalous state occurs. A complete state diagram includes 
both normal and anomalous states. The transition probability matrix 
is implicitly given in Figs. 2(3) and 4 for rate \ (§) codes. The Markov 
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Chain is regular and therefore the stationary probabilities exist. Let 
Pi , i = 1, 2, ■ • • , 2 0+l , A , Ai , • • • , A v be the stationary probabilities 
of being in the respective states. 

An upper bound on the probability of block error, P t „i , is obtained 
by summing the stationary probabilities that the decoder is in a normal 
state whose last digit is 1 or in an anomalous state A through A v -i • 
If the decoder reaches state A v then correct decodings are assured in 
states A A + 1 to A v , but if the decoder returns to state A after state 
A A + i , j = 1, 2, • • • , (V — A — 1) correct decodings are not assured 
in states A A + i to A K + j . A tighter bound which accounts for this dif- 
ference will be obtained next. For rate f codes with t = 1 and d m ^ 4 

Peul — P5 + Pfi + P8 + Po + PlO 4" Pl3 + Pl4 + Pl6 + PA. 

+Pa> + Pa, + Pa, + Pa. • (13) 

Expressions for F, ul for rate \ and f codes are given in Table I. 

A tighter upper bound, P eu2 , is obtained by using 2V anomalous 
states as in Fig. 5. State B exists whenever channel errors occur. 
States Bi , B 2 , • • • , B v -i represent successive blocks free of channel 
errors when the total number of consecutive blocks free of channel 
errors is less than V. States C, , C 2 , ••• ,C V represent successive blocks 
free of channel errors when the total number of consecutive blocks 

Table III— A, A, V For Typical Values of Code Parameters 
Code Parameters A A 1 
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1-Po 




Pd(Po+Pw) bH 



Fig. 4 — (V + 1) anomalous states for P mi bound on probability of block error. 



Po^2 Pofo Po/?4 




Pd(Po+Pw) b "' 



/?L = l-Po V " 1 1 = 1,2,... 

Fig. 5 — (2V) anomalous states for P eu2 bound on probability of block error. 
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free of channel errors is at least V. In Fig. 5 state X represents any 
normal state marked X in Figs. 2 and 3 from which the transition 
from normal state to anomalous state occurs. A complete state diagram 
includes both normal and anomalous states. The transition probability 
matrix is implicitly given in Figs. 2(3) and 5 for rate ^(f) codes. The 
Markov Chain is regular and therefore the stationary probabilities 
exist. 5 Let s { , i = 1,2, • • • , 2 + , B , B, , • • • , B v -- L , Ci , C 2 , • • • , C v 
be the stationary probabilities of being in the respective states. 

P eu2 is obtained by summing the stationary probabilities that the 
decoder is in a normal state whose last digit is 1 or in an anomalous 
state B through B v -\ or Ci through C A . For rate f codes with t = 1 
and d m ^ 4 

P.u2 = S 5 + S + S s + S + S 10 + S 13 + Sm + S, 8 -f- S„ 

+ S fll + Sb, + s Bl + s Bt + s Cl + s c , . (14) 

Expressions for P„ u2 for rate \ and f codes are given in Table I. 

5.2 Lower Bound on the Probability of Block Error 

A lower bound on the probability of block error, P, L , is obtained 
by summing the stationary probabilities that the decoder is in a normal 
state whose last digit is 1. That is, it is assumed that no decoding 
errors are made when the decoder is in anomalous states. For rate 

§ codes with t = 1 and d m ^ 4 

P.L = PS + P0 + P8 + P9 + PlO + Pl3 + Pli + Pl6 • (15) 

Expressions for P eL are given in Table I. 

VI. PERFORMANCE OF TWO CODES 

The performance of a rate \ (18, 9) code of minimum distance d m = 6 
has been evaluated on the binary symmetric channel {I = 1). This 
code is obtained by extending the (17, 9) quadratic residue code of 
minimum distance d m = 5 by one bit. Upper and lower bounds on 
the probability of block error for the (18, 9) code used in burst-trapping 
with t = 1 and t = 2 versus the binary symmetric channel transition 
probability, p, are given in Fig. 6. The probability of block error for 
maximum likelihood decoding of the (18, 9) code and the probability 
of block error for an uncoded 9 bit block are presented for comparison. 
It is interesting to note that the performance of the code is better 
with burst-trapping decoding (t = 2) than with maximum likelihood 
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BLOCK ERROR RATE 



Fig. 6 — Performance of (18, 9) code on binary symmetric channel, a. Upper and 
lower bounds for burst-trapping procedure, t = 2; b. Upper bound for burst-trapping 
procedure, t ■» 1; c. Maximum likelihood decoding; d. Uncoded 9-bit block; e. Deci- 
sion feedback, t = ; f . Decision feedback, t = 2. 



decoding for values of p where the rate 5 is much less than the channe 
capacity. Channel capacity is | when p = 0.11. The probability of 
block error for two decision feedback systems using the (18, 9) code 
are also presented for comparison. In one system, (t = 0), repeat 
requests are made on all received n-tuples but code words. In the 
other system, (t = 2), single and double error corrections are made 
and repeat requests are made on all received n-tuples but code words 
or code words perturbed by one or two errors. 

The performance of a rate \ (39, 26) shortened BCH code of minimum 
distance d m = 6 with t = 1 has been computed on the basis of recorded 
telephone error data. One set of recorded telephone error data is the 
Vestigial-Sideband (VSB) data. 6,7 A selected set of 85 calls was used. 
Each call is an error sequence of about 3 X 10 6 bits recorded at 3600 b/s 
(4-level operation of the VSB modem). The 85 calls were divided into 
groups of calls of similar bit error rate as shown in Table IV and P(m, 39) 
statistics 8 were determined for each group. P(m, 39), m = 1,2, • • • ,39, 
is the probability that m errors occur in a block of 39 bits. Estimates 
of p , p e , p K , p d , p u , and a were obtained for each group by using 
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approximate code and coset weight spectra. A plot of the distribution 
of the percent of calls (quantized by the grouping of calls) with respect 
to raw (39 bit) block error rate and the upper and lower bounds on the 
distribution with respect to decoded block error rate is given in Fig. 7. 
The distribution of the percent of calls with respect to decoded block 
error rate as estimated by computer simulation 2 on the same 85 call 
sample is also given for comparison. For the comparison simulation a 
block interleaving degree t = 117 was used as this amount of inter- 
leaving attained approximate independence between successive blocks 
of the subcodes. 

Another set of recorded telephone error data is the Alexander-Gryb- 
Nast (AGN) data. AGN data consists of about 1000 calls recorded 
at data rates of 600 b/s and 1200 b/s on an FM data modem. These 
calls are from three classes of calls of approximately 3.6 X 10 5 , 7.2 X 10 6 , 
or 2.16 X 10° bits in length. The 294 calls with bit error rate greater 
than 10~ 5 were divided into groups as shown in Table V and estimates 
°f Po , Pc , P« t Pd } p u , and a were obtained as for the VSB data. A 




BLOCK ERROR RATE 



Fig. 7 — Distribution of 85 call VSB sample with respect to block error rate. a. 
Uncoded 39 bit blocks; b. Simulation of (39, 26) code interleaved to degree 117; 

c. Upper and lower bounds for (39, 26) code. 
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plot of the distribution of the percent of calls (quantized by the group- 
ing of calls) with respect to raw (39 bit) block error rate and the upper 
and lower bounds on the distribution with respect to decoded block 
error rate is given in Figs. 8 and 9. The performance as estimated by 
computer simulation 2 on the same call sample is also presented for 
comparison. 

For both sets of data there is good agreement between the performance 
estimates obtained by the theory and simulation. Two advantages of 
the theoretical technique over the simulation technique are the ability 
to obtain the tail of the distribution and the ability to obtain the 
distribution from reduced data. This latter advantage is due to the 
fact that p , p c , Pw , Pd , p u , and a can be approximately determined 
from reduced data such as P(m, n) statistics. 

VII. CONCLUSIONS 

A technique for analyzing the performance of a burst-trapping error 
correction procedure has been described. The criterion of performance 
is the probability of block error. The analysis technique for burst- 
trapping procedures with block interleaving degree I is valid on channels 
where error patterns in blocks separated by t blocks are (approximately) 
independent. Thus the analysis technique is valid for random error 
channels or for properly designed burst-trapping procedures on burst 
or compound channels such as the telephone channel. 

The performance of a rate f code, the (39, 26) shortened BCH code, 
has been computed for two sets of recorded error data. For both sets 
of data the computed performance agrees with the performance obtained 
by computer simulation. An advantage of the theoretical technique over 
simulation is the ability to estimate performance on the basis of reduced 
data [P(m, n) statistics] rather than on extensive error sequence data. 
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APPENDIX A 



Transition Probabilities for Normal State Diagrams 

For simplicity the development of the one step transition probabilities 
p iti from renumbered state i to renumbered state j is demonstrated for 
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Fig. 8 — Distribution of AGN call sample with respect to block error rate. a. 
Uncoded 39 bit blocks; 2. Simulation of (39, 26) code interleaved to degree 117; 
c. Upper and lower bounds for (39, 26) code. 

rate f codes (Fig. 3). The technique is easily generalized to rate (b— \)/b 
codes, 6^2. Obviously, 



Pi.i = Po + p c ; 
Pi. 9 = p u + p» ; 

Pi. 2 + Pi. 10 = P,l ■ 



(16) 
(17) 
(18) 



Let Pi,i(b) be the b step transition probability from renumbered state 
i to renumbered state j. Due to the structure of the normal state diagram 
it is easily seen that, 



Pi.i = P.. 4(3) + p,. 5 (3) + p ll7 (3) + p ll8 (3); 

,, v Pi.4(3)+Pi, 5 (3) . 
P,., = (1 - p 2 .,) = , 

Pi .2 

n \ Pi ■< (3) . 

*■< = ° " P3 - s) = pT^; • 



(19) 
(20) 

(21) 
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10 ° 10"* 

BLOCK ERROR RATE 



Fig. 9 — Distribution of AGN call sample with respect to block error rate. a. 
Uncoded 39 bit block; b. Upper and lower bounds for (39, 26) code. 



Pa. v = (1 — Pe.s) = 



Pl.7(3) 

Pi. 2^2,6 



(22) 



The three step transition probabilities pi. 4 (3), Pi. 5 (3), Pi. 7 (3), Pi,s(3) 
can be computed by summing the probabilities of all combinations of 
three successive error patterns which cause transition from normal 
state 10 to normal states 4000, 4001, 4010, and 4011 respectively. 



Pi, 4 (3) = p d p 2 o ; 

Pi. s(3) = VdVoPv ; 

Pi. ?(3) = p rf p tt po ; 

Pi. 8 (3) = VdVl • 



(23) 
(24) 

(25) 
(26) 



Further, 



Pi.io = P,., 2 (3) + Pi., 3 (3) + Pi., 5 (3) + Pi,, 6 (3); 



(27) 



PERFORMANCE OF BURST-TRAPPING CODES 513 

Pio.u = (1 - Pio.u) - - f W 

Pi. 10 

p„... = (i - p.,„) = -f^ ; (29) 

f'l.lO/'lO.ll 

P„... = (1 - P„.,«) = ^ 5L - (30) 

fl, 10/^10,14 

The three step transition probabilities pi.i 9 (3), Pi.i 3 (3), Pi,is(3), Pi.ig(3) 
can be computed by summing the probabilities of all combinations 
of three successive error patterns which cause transition from normal 
state 10 to normal states 4100, 4101, 4110, 4111 respectively. 

P.. i 2 (3) = p d {a(p c + p u + pS~ + 2p a(p c + p u + p d )\ ; (31) 

Pi. i3(3) = p.,{a(l - a)(p e + p„ + p d f + Po(l - a)(p e + p u + pa) 

+ a(p c + p u + Pd)p«.} ; ( 32 ) 

P,.i 5 (3) =Pi, 13 (3); (33) 

Pi.mO) = P-t(l - «) 2 (Pc + P« + P^) 2 + 2 P-(! - «)(?>< + P- + P*)l • 

(34) 

The transition probabilities for rate § codes are given in Fig. 2. The 
transition probabilities for rate f codes with the exception of 

a(Pc + Pu + p,i) + 2ap , ,o,r\ 

*»•» = «^ ' c ; 

(1 — a) a(p c + p u + Prf) + Po(l — a) + «P«, . f oa\ 

P "-" = « + Po ' ^ 

(1 — a) a(Pc + Pu + Prf) + Po(l — g) + op.- . /Q 7 \ 

Pl4 ' 15 = (1 - «) + p. ' (d?) 

(1 - a)\p c + p„ + P„) + 2(1 - a)p w . . v 

Pl4 ' 16 = (1 - a) + P. ' (38) 

are given in Fig. 3. 

APPENDIX B 

Bounds on the Amount and Length of Error Propagation and Determina- 
tion 0] the Recovery Space 



514 THE BELL SYSTEM TECHNICAL JOURNAL, APRIL 1970 

We consider codes in which d m ^ 2t + 2. Let ^ e t ^ k be the 
maximum number of errors possible in the information part of the 
decoded word at block [i€[. The value of e, depends on the method 
of decoding the block and the error pattern sequence. We will say that 

ft = Z e,- (39) 

is the maximum number of errors "available" for feedback to blocks 
after block [it]. As in the text, assume that one or more of blocks [(i — 
b + 2)1] to [U] is incorrectly decoded and that block [ii] has channel 
errors, but all succeeding blocks [(i + l)f\, [(i + 2)1], • • • , have no 
channel errors. We consider three mutually exclusive and exhaustive 
cases. We use the fact that a block code of minimum distance d m is 
simultaneously capable of correcting t errors and detecting d ^ t 
errors if d m = t + d + 1. 

b.i Case 1 

Assume that no error detections (BTD's) occur in blocks [(i + l)t], 
[(i + 2)t], • • • . Under this hypothesis error propagation cannot occur 
beyond block [j(] if /,• < d m — t. The maximum length of error propaga- 
tion is then the smallest integer Ai such that /j+Ar < d m — t. Since 
/,• = (b — l)k, blocks [(i + 1)£] to [(i + b — l)t] are subject to decoding 
error due to error propagation and therefore /,+6-i = (b — l)t. Each 
subsequent block, [{i + b)(\, [{i + b + 1)1], • • • , which is correctly 
(incorrectly) decoded removes t (d m — 2t) errors from the maximum 
number available for feedback. Therefore 

'(*-« + If^jl. ifd.-*S< 

A, = J Ld ~ ~ 1U (40) 



(b _ 1) + | M - 2Q. - 20 j | 



if d m - 2t > I 



where [x_\ is the greatest integer less than or equal to x. 
The maximum amount of error propagation is 

(b - 2)t 



A 1 = (6 - 1) + 



d m - 2t 



(41) 



Blocks [{i + Ai + b)l], [(i + A, + b + 1)1], • • • are correctly decoded 
by RED in normal state 10. 

b.2 Case 2 

Assume that an error detection (BTD) occurs in block [(i + j)t], 
1 ^ j == b — 1. Under this hypothesis error propagation cannot occur 
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beyond block [(i 4- j)4- The maximum length of error propagation is 

A 2 = 6 - 1 (42) 

and the maximum amount of error propagation is 

A 2 = b - 1. (43) 

Blocks [(i + A 2 + b)4, [(* + A 2 + 6 + 1)4, • ■ ■ , are correctly de- 
coded by RED in normal state 10. 

b.3 Case 3 

Assume that an error detection (BTD) occurs in block [(i + j)t), 
b ^ j ^ J, where J is the maximum value of j such that block [(i 4- i)4 
can be detected in error when the last channel errors occur in block 
[U] and blocks [(* + 1)4, [(* + 2)4, • • • are free of channel errors. 
Under this hypothesis, error propagation cannot occur beyond block 
[(* + i)4- For block [(* + J)(\ to be detected in error it must have 
k ^ * + 1 errors fed back to it. For 6 = 2 (rate § codes) J = b — 1, 
and Case 3 cannot occur. 

As in Case 1, /,- = (6 - l)fc, and blocks [(t + 1)4 to [(t + 6-1)4 
are subject to decoding error due to error propagation. Therefore, 
/,+»_, = (6 - 1)*. Each subsequent block, [(» + 6)4, [(* + & + 1)4, 
which is correctly (incorrectly) decoded removes <(d m — 20 errors from 
the maximum number available for feedback. Let H a be the largest 
integer such that f i+Ha ^ t + 1. Then when 6^3 



tf „ = 



(&-!) + 



(6 — 2)t — 1 

d m - 2< 



, if d n - 2i ^ i; 

(44) 
26-4 if d m - 2t> t. 



Since errors must be fedback to block [(» + J) 4 from at least two 
blocks to get t x > t feedback errors, 

,/ = H a 4- (6 - 2). (45) 

b.3.1 Case 3a 

Assume block [(i + ;*)4 is correctly decoded by BTD. Then 



A 3 = H a 






(46) 


A 3 = (6 - 


- D + 


1 (6 - 2)1 - 1 

L d m - 2t _ 


(47) 



Blocks [(* + J + 6)4, [(i + J + 6 + 1)4, * ' ■ are correctly decoded in 
normal state 10. 
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B.3.2 Case 8b 

Assume block [(i + j)l] is incorrectly decoded by BTD. Now b ^ 
j ^ G ^ J, where G is the maximum value of j such that block [(i + j)l] 
can be detected in error and incorrectly decoded due to error propaga- 
tion when the last channel errors occur at block [U] and blocks [(i + l)t], 
[(i + 2)1], • • • are free of channel errors. For block [(i + G)(] to be 
detected in error it must have t x ^ t -\- 1 errors fed back to it to cause 
the detection and * 2 = 1 errors fed back during the BTD to cause er- 
roneous decoding. 

If /i + j,_i = (6 — l)t < t + 2, (7 = b — 1 and Case 3b cannot occur. 
Let H b be the largest integer such that /,- +J / t > £ + 1. Then when 
6 ^ 3 and (6 - 2)t ^ 2 



#„ = 



L rf m - ^ J (4g) 

(6 _ 1) + I fl ~ 2 >* ~ 2 1 , if rf m - 2i > t. 



Since errors must be fed back to block [(i + G)l] from at least two 
blocks (three blocks if t = 1) to get t t + t 2 ^ t + 2 feedback errors 

A4 = G ={ ff ' + (»- 2 >- if '*'! (49) 

Ifl, + (6 - 3), if ( - 1. 

^..p-D+ ^-q' -« + '--* ]. (so. 

Blocks [(i + A 4 + &)/], [(z + A 4 + & + l)fl, ■ ■ ■ are correctly decoded 
in normal state 10. 

Table I is constructed by taking 

A = max {A x , A 2 , A 3 , A 4 } (51) 

and 

A = max [A t , A 2 , A 3 , A A ). 

By definition the recovery space is the minimum number of con- 
secutive blocks free of channel errors required to guarantee the return 
of the decoder to one of normal states { 10, 11, 20, 21} at block [{i + 
y+1)^] regardless of the channel error sequence prior to block [(i+l)fl. 
We say that V is a necessary recovery space if fewer than V consecutive 
blocks free of channel errors do not guarantee transition to one of 
normal states J 10, 11, 20, 21} at block [(t + V + 1)1]. We say that 
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7 is a sufficient recovery space if V consecutive blocks free of channel 
errors guarantee transition to one of normal states (10, 11, 20, 21} 
at block [(i + V + l)t\. For Cases 1, 2, and 4 7,- = A, + (b - 1), 
i = 1,2, and 4, is a necessary and sufficient recovery space. For Case 3 
7 3 = A 3 + (26 — 3) is a necessary and sufficient recovery space. Since 
in all cases V ^ V, , i = 1, 2, 3, 4, is a sufficient recovery space V = 
max [Vi , V 2 , V z , V*) is a necessary and sufficient recovery space. 

APPENDIX C 

List oj Symbols 

I — block interleaving constant 

b —parameter determining rate of recurrent code 

A — recurrent code parity check matrix 

H — component block code parity check matrix 

N — constraint length of recurrent code 

A N —recurrent code truncated parity check matrix 

n — block length of component block code 

k — number of information bits per block of component block code 

I it — /c-tuple of information bits of [U]t\\ block 

M it — n-tuple of transmitted bits of [U]th block 

Qit — ( n — fc)-tuple of parity bits of [i(\th. block 

I\, — jth segment of In 

M& — n-tuple of received bits of [il]t\\ block 

I* —/c-tuple of received information bits of [U]th block 

Q* — (n — fc)-tuple of received parity bits of [U]th block 

P*c — parity (n — fc)-tuple obtained by encoding I? t with H 

I it — /c-tuple of decoded information bits of [it\th block 

P'ic — Q* as modified by fed back information segments 

I' it — jth segment of In 

Co — the set of one element— the pattern with no errors 

C w — the set of channel error patterns which are identical to nonzero 

code words 
C c — the set of nonzero channel error patterns which are correctable 

by RED 
C„ —the set of channel error patterns which are both uncorrectable 

by RED and undetectable 
—the set of channel error patterns which are detectable 
— the probability that the information bits of a received word 

are unaffected by the channel error pattern given that the 

channel error pattern is correctable, detectable, or undetectable 



a 
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p„ — the probability that the channel error pattern is in set C 

p w — the probability that the channel error pattern is in set C w 

-p e —the probability that the channel error pattern is in set C c 

p u — 'the probability that the channel error pattern is in set C u 

■p d — the probability that the channel error pattern is in set C d 

d m — minimum distance of the component block code 

t — amount of error correction done in RED 

V ■ — -recovery space 

p,,, — single step transition probability from state i to state ; 

7T, — stationary probability of being in state i with Genie Decoder 

P, tt — probability of block error with Genie Decoder 

a ■ — amount of error propagation 

A — upper bound on a 

X — length of error propagation 

A — upper bound on X 

A, — designation of anomalous states for upper bound one on block 

error probability 
Pi — stationary probability of being in state i for upper bound one 
Peu\ — upper bound one on probability of block error 
B { ,d — designation of anomalous states for upper bound two on block 

error probability 
Si ■ — stationary probability of being in state i for upper bound two 
P*„2 — upper bound two on probability of block error 
P tL ■ — lower bound on probability of block error 
Pi. i (3)— three step transition probability from state i to state j 
e, — maximum number of errors possible in decoded information 

fc-tuple at block [U] 
fi — maximum number of errors available for feedback to blocks 

after block [it] 
d —amount of error detection capability of component block code 
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